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Abstract 



Heat conduction in ionized plasmas in the presence of magnetic fields 
is today a fashionable problem. The kinetic theory of plasmas, in the con- 
text of non-equilibrium thermodynamics, predicts a Hall-effect-like heat 
flow due to the presence of a magnetic field in ionized gases. This cross ef- 
fect, the Righi-Leduc effect together with a heat flow perpendicular to the 
magnetic field are shown to yield significant contributions, under certain 
conditions to the ordinary Fourier component of the heat flow. The ther- 
mal conductivities associated with these effects change with the strength 
of the magnetic field for a given temperature and density and are shown to 
be significant compared with the parallel conductivity for a whole range 
of values of B. 



1 Introduction 

The behavior of charged particles in the presence of external magnetic fields 
is a well known subject. According to the tenets of classical non-equilibrium 
thermodynamics for not too large gradients, the simultaneous heat and electrical 
flows in a magnetic field are linear functions of the temperature and electrical 
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potential gradients, respectively. The coefficients appearing in these relations 
are, in general, second rank tensors which are functions of the external field 
and consist of a symmetrical part and an antisymmetric Al part [HOIS]. The 
latter ones are called the "Hall vector" in the case of electrical conduction and 
the "Righi-Leduc" vector in the case of heat conduction. The former is a well 
known effect whereas the latter one isn't. It was discovered in 1887 and first 
measured and confirmed by Waldemar Voigt in 1903 (see Refs. [HE]). Curiously 
enough it has been, ever since, hardly mentioned in the literature. Not even 
the authors of a rather beautiful and striking experiment published recently |4J 
recognize that what they have really detected is the Righi-Leduc effect. 

In this paper we address ourselves to a rather different calculation of the heat 
conduction effects in magnetized plasmas from the theoretical point of view. 
Indeed, in the case of a dilute ionized gas in the presence of weak magnetic 
fields for densities n in the interval 10 3 < n < 10 s m~ 3 and temperatures in the 
range 10 3 < T < 10 7 K the magnitude of heat, charge and mass currents arising 
from these and other so-called "cross effects" may contribute by factors which are 
not necessarily negligible when compared with those arising from ordinary heat, 
mass and electrical conduction. This result seems to differ from those obtained 
by Braginski [5] while qualitatively agree with those obtained by Balescu in his 
excellent and exhaustive treatment on plasma transport processes [6]. We shall 
come back to this point a little later in this paper. Moreover, we believe that the 
approach here offered to the general subject of collisional transport processes 
in plasmas may provide a better understanding of this subject at temperatures 
and densities where collisions are favored in a non-relativistic framework. 

The structure of the paper is as follows. In section II we shall briefly review 
the kinetic concepts used throughout the work, section III contains the bulk of 
the main calculations leading to the heat conduction transport coefficients and 
in section IV the results are given together with some concluding remarks. 

2 Kinetic model 

The kinetic model we use in our approach to the problem is the well known 
one based on the Boltzmann equation as originally proposed by Chapman and 
Cowling [7]. The substantial difference is that for the case of ionized gases they 
never developed their method up to the stage of placing it within the framework 
of linear irreversible thermodynamics. Although they derive elementary expres- 
sions for the Righi-Leduc and the Nernst-Ettingshausen effects (see Ref. [2]) 
they never pursued any comparison of their formulae with experiments, much 
less studied any possible applications. Here we would like to focus ourselves on 
the former and related effects. 

Thus, the system here considered is a binary mixture of electrically charged 
particles with masses and charges m, and e% for i — a,b. The density of the 
system is low enough such that the kinetic description is valid. For simplicity 
we shall set the charge of the ions Z = 1. For such a system, the evolution of 
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the distribution functions of the molecules is given by the Boltzmann equation: 



where the subscript i indicates the species and J (fifj) is the collisional term 
representing collisions between different and same species. The Lorentz force 
on the left side of Eq. ([1]) is treated here as an external force where the mag- 
nitude of the magnetic field is small enough such that collisions dominate over 
cyclotron motion. We must clarify that both the electric and magnetic fields 
appearing in this force contain the self consistent fields produced by the plasma 
and are governed by the Maxwell's equations (see Ref. [6]). This weak field 
approximation implies u>iT ps 1 where u>i — eB/rrii is the Larmor frequency, 
that is, the frequency of the circular orbits that the particles describe around 
magnetic field lines. 

Once Eq. (I]) is defined, the derivation of the conservation equations as well 
as the proof of the H-theorem are logically required. The former is a standard 
step widely discussed in the literature (|6j-[10j). The second result becomes 
a little bit more tricky due to the presence of a magnetic field but we will 
not bother with it since it is not essential to this paper [8]. We here proceed 
directly with the solution of Eq. (T} following the standard Hilbert-Chapman- 
Enskog approximation. Since a local Maxwellian distribution function is 
clearly a solution to the homogeneous part of Eq. (U) we propose that the single 
particle distribution functions fa (i = a, b) may be expanded around in a 
power series of Knudsen's parameter e which, as well known, is a measure of the 
magnitude of the macroscopic gradients [7l [9J, [TO] . Thus 



In Eq. ((2j) the functional equilibrium assumption is also invoked, namely the 
time dependence of f+ and ip^ for all n occurs only through the conserved 
densities. The particle density is rii(r,t) (i = a, b), the barycentric velocity 
u(f 7 t) and the local temperature T (r,t) is in this paper assumed to be the 
same for ions and electrons [IT]. T (r, t) is related to the internal energy density 
e (r, t) by the standard ideal gas relationship. In this work we shall deal only 
with the first order in the gradients correction to fi characterized by tp^ , namely 
the Navier-Stokes-Fourier regime. 

Substitution of Eq. <[2j) into Eq. (TJ), leads to order zero in e to the Euler 
equations of magnetohydrodynamics. To first order in e one obtains a set of two 
linear integral equations for ip^ which involve the linearized collision kernels 
in their homogeneous terms whereas the inhomogeneous ones contain a combi- 
nation of terms involving the macroscopic gradients VT, Vu and the diffusive 
force dij — —dji. Indeed, after a somewhat lengthy but standard manipulation 




(1) 




(2) 
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[7] one gets that 



m^c 2 5\ VT Uj T n 
a 



H^j 

2kT 2] T n" u 



i = a, b 



^ x + c (rf>) + c(,fVf) 



(3) 



where a superscript "o" over a tensor indicates its symmetric traceless part. In 
Eq. J3]), C (vi^) an d C f^Vf^J are the linearized collision kernels whose 

explicit forms are also well known [6]-[l0] and dij is the diffusive vector force 
given by 

7 , n a n b (m a - m b )Vp n a n b -> 

d a b = V 1 ym h e a - m a e b ) ■ E' (4) 

n np p pp 

which satisfies the property d^ = —dji. E' = E + u x B is the "effective" electric 
force, u the bary centric force defined as 



pu 



(r,t) = J2pA{r,t) (5) 



Pi = mini (i = a, b) and p = p a + p b . Also u, (r,t) — (3) where () is the 
standard average taken with fi (r, Vi,t). 

We must emphasize that, as has been recently shown [10J, it is only in this 
representation that the Onsager reciprocity relations hold true, at least when 
B = 0. If B is different from zero the proof of such relations starting from the 
linearized integral equations for tp^ will be discussed elsewhere. The solution 
to Eqs. (J3j) has been carefully outlined in Ref. [13]. Using Curie's theorem 
which allows setting Vu = 0, the solutions are found to be of the form 

= ■ -jT - &i ■ dy i = a,b (6) 

where 

A t = A^Ci + A< 2) (ci x B) + Af ] B (c* • b) (7) 

and a similar expression for Dj. The scalars Ap } (j = 1, 2, 3) appearing in Eq. 

(7]) are functions of of n, T, cf , -B 2 and ^3 ' ^ • 

When Eq. |(6|) and its analog for D,. are substituted in Eq. ([3]) we get a set 
of linear integral equations for the unknown functions A^ and their analogs for 
the Dj vector. As shown in Ref. (13], these equations are 

if (|gl - 5) 3 = /i 0) {c (3^) + c {mi + 3^)} (8) 
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where IZi = A 



_ a(i) 



B 2 A 



(3) .•_ 
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p(0) 
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2fcT 
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In Eq. 
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iBclG-f 
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(9) 
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(i) 



iBA 



(2) 
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(i) 



iBG 



(2) 



1 b 



<^/j 0) Af 



1 

B 2 



(10) 
(11) 

(12) 



for k = 1, 2. 



(7) 

Similar results are obtained for the unknown functions D ■ ; appearing in the 
analog of Eq. (7} for the vector but we shall not bother with them here since 
we will concentrate only in heat conduction in the plasma. The ^vector is 
clearly related to diffusion phenomena [8]. Equations and © are the basic 
ingredients required to discuss the problem of heat conduction. The curious 
reader will immediately realize the difference between these basic results and 
those used by Braginski in his paper on this subject [Ref. [5] pages 245-247]. 
In our method the unknown Aj functions satisfy integral equations, ((8| and 
iJH), where the collisional dynamics obeyed by the particles of the same and 
different species are contained in the linearized collision kernels about which 
no assumptions have yet been introduced. We shall now proceed in the next 
section to study the heat conduction in the dilute plasma. 



3 Heat conduction in a fully ionized plasma 

According to classical irreversible thermodynamics [TJ EE [14] the expression for 
the heat flux vector in a multicomponent system is given by 

Jq = Jq ~ / . hi 
H rrn 

1 

where Ji is the diffusion vector for species i and hi its enthalpy. Since for an 
ideal gas hi = \kT where k is the Boltzmann constant, using the standard 
definition for J q , namely for our case 



i— a 



we get that 

kT" q ~ -fW V 2kT 2 
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since Ji = mi 



Also, J' = for fi — ff^ so that, substituting Eq. ([6]) in 



Eq. lfT3|) and ignoring diffusive contributions to J' q we get that 



(fcT)" 



4 = E 



2fcT 2j /l ^ l Cl<Z l 



(14) 



where ^j 1 ' is now given by Eqs. and J7]). 

To continue with the calculation of J' q we now resort to the almost orthodox 

* (k) 

method in kinetic theory, namely to expand the unknown functions Aj in 
terms of a complete set of orthonormal functions, the Sonine polynomials, so 
that 



E«. 

rn—0 



(fc)(m) q(m) i 2\ 
3 °3/2 \ C i) 



k = 1, 2, 3 



(15) 



where the coefficients ^ are still to be determined from the integral equa- 
tions (HI and ([9]) and thus depend on the interaction (Coulomb) potential be- 
tween the species. Notice however that in spite of the complicated form for 
Aj given in Eq. ([7]), when Eqs. <(6j) and lO are introduced into Eq. (fl"4"l) all 
integrals over Ci have the same structure, namely 



E 



2kT 



(0) a (k)(r, 
m=0 



'^3/2^ ,( 



where use has been made of the well known property that |6J-[10J 

„2n+l dx _ r(n+p + g ) f 



and that Sn {x) 



~* S« (x)5^ Or)* 2 
1. Clearly then 



2p 



fc= 1, 2, 3 
(16) 

(17) 
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rri; 



« (1)(1) VT- 



of» (1 »Bx VT- 



. (3)(1) B(S . y T 



(18) 



Equation (|18|) is an important result. Of the infinite number of coefficients 
required to specify the functions A^ appearing in Eq. ([15)) only one is required 
to compute the explicit form of the constitutive equation <fl8|) . This fact readily 
simplifies the solutions to the integral equations ([8]) and |(9|) as we shall see 
below. However, before doing so let us examine Eq. (fl8|) . If we take B in the 
direction of the z-axis, B x VT is a vector perpendicular to both B and VT 
whereas B (j3 ■ VT^j = B 2 ^k so that we may write that 



b 

E 



k 2 T^^ 



(af )(1) 



i?V (3)(1) 



ViiT- 



a, (1)(1) V,T 



Baf ){1) V s T 



(19) 



6 



where the last term represents a heat flow in the direction perpendicular to both 
B and VT. This is the well known Righi-Leduc effect [I]-[3]. Notice also that 
when B = 

4— |^E£M 1)(1) VT = - K|| VT (20) 



i— a 



which is the well known form for Fourier's heat conduction equation and the 
thermal conductivity is given by 



«ll)B=o = !^E^f )(1) (21) 



i— a 



In Eq. l[2~Tj) the coefficients ai ^ and a^ 1 must arise from the solution to 
the integral equations and which indeed become identical to each other 
when B — 0. In its more general form, Eq. (fT9| thus reads as 

f q = -«||V||T - K^VxT - k s V s T (22) 

where the three conductivities are readily identified from Eq. lfl9|) . Once more 
a word of caution. Equation lj22|) appears to be the same as the one quoted by 
Braginski [see Eq. (4.33) of Ref. [5]] however, the evaluation of the coefficients 
a (!)( m ) m our cage rac jically differs from the procedure followed by this author. 
And worst, his definition of heat flux is foreign to the one used here, that is, not 
in agreement with irreversible thermodynamics. 

To determine the ^ coefficients we need to solve Eqs. and ((9]). 
The former one is straightforward and has been discussed in the literature, 
specially in Appendix B of Ref. [13]. The second one, however poses some 
problems due to the structure of its inhomogeneous term. Nevertheless by a 
subtle generalization of the procedure followed to solve Eq. ([8]) , Eq. J9]) is also 
solved using a variational procedure apparently due to Davison [l5j[T6]. The 
steps are also outlined in the same Appendix B of Ref. [13| . We shall only quote 
the results here, 

a^+B 2 a^ =2.94r 

a««+B 2 ai 3)(1) = l-96r (23) 

Ji)(i) , R2„(3)(l) _ °- 058 _ 
% +Ba b -—r 

also calling af 1 ^ = a^^ m ^ + iBaf^ m \ i = a, b, m = 1, 2 following Eqs. IjlOp 
and (fl5| we get, after separating real and imaginary parts that 



Re 



ai 0) 



aW (0) = (56.3 - 662x 2 - 2.25x 4 ) 



Im 



= Ba^ = (647x + 2.2x 3 ) 
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Re 
Im 



= ai 1 ^ = (37.5 + 2U7x 2 + 7.3a; 4 ) 4~ 
= Bai 1){2) = (206a; + 2649a; 3 + 9a; 5 



T 

A, 

T 

a7 



(24) 



Re 

Illl 



= (1.12 + 121. 2x 2 + 154.4x 4 



(1)(2) 



(0.06a; + 7a; 3 + 9a; 5 



MiAi 

T 



where 



Ai = 19 + 2078a; 2 + 2650a; 4 + 9a; 6 



(25) 

and x = ui e T = 1.76 x 10 11 Bt where B is given in teslas. r is the mean free 
time obtained from the only independent collision integral and is defined a£] 



4(27r) 3/2 y77i;(fcr) 3/2 eg 



(26) 



In Eq. (|26|) ip is the so called logarithmic function which arises from using 
Debye's length as a cutoff length in the collision integrals, defined as 



ip = In 



1 



/ 16tt£;TA 



where Xd is the Debye's length, namely 



A, 



e Q kT 



(27) 



(28) 



B=0 



= (Re[a! 13 ]) 



B=0 



Notice that the consistency requirement that ^ai 1 ^ 1 ^ 

satisfied within the approximation used here. 
Therefore, summarizing the three thermal conductivities in Eq. (|22|) are given 
by 

5 k 2 T 

K\\ = x 2.01tit (29) 



4 m R 



; = ^ — ^ x (38.7 + 2270a; 2 + 161a; 4 ) 



4 m e Ai 
5 k 2 T nr 



(206a; + 2644a; 3 



(30) 
(31) 



4 m e Ai 

where ip is defined in Eq. (|27|l . 

We emphasize that Eqs. (f29|) - (f3T|l are valid for a fully ionized hydrogen 
plasma (n a — rib — n/2, mi, 3> m a — m e ) in a first order in the gradients 
approximation, namely the Navier-Stokes-Fourier regime. These are the main 
results of this paper. 



lr The numerical factors appearing in the six polynomials quoted in Eqs. I|24|l substitute 
those published in Appendix B of Ref. |13| which result from a revised calculation. 



s 



4 Discussion of the results 



The first thought that a reader may have concerning the nature of Eqs. (|29"ll - 
lf3lj) is to ask how they compare to those obtained from the well accepted and 
dominant calculations obtained much earlier by Spitzer [17] and by Braginski [5] . 
As already pointed out in full detail by Balescu in Ref. [6] and emphasized by us 
in a recent paper [T3| this is quite difficult. Neither of both authors, besides using 
the Fokker-Planck and Landau kinetic equations, respectively, performed their 
calculations within the framework of classical irreversible thermodynamics. In 
particular, the correct definition for the heat flux in a multicomponent mixture, 
e. g. Eq. (fl2|) . was ignored. Therefore the expressions that they quote for the 
thermal conductivities are not equivalent to ours. 

The most interesting feature of our results is that the three thermal conduc- 
tivities, exhibit a behavior which is similar to that shown by the results obtained 
by Balescu who used a Landau type kinetic equation which he solved using the 
"moments" method. This result is gratifying. It shows indeed that using the full 
Boltzmann equation and solving the ensuing integral equations which define the 
coefficients appearing in the explicit results for the transport coefficients leads 
to results that approximately equal to those obtained from the Landau equa- 
tion when the practically "exact" 21 moment approximation is used. Comparing 
figure 1 of our paper with figure 5.1 in Balescu's book confirms this statement. 
Thus, his prediction that different methods used to compute transport coeffi- 
cients will yield results that will be within a 10% difference from each other 
turns out to be sustained. 

Another issue is important. Working with the moment method leads to 
results which are not clearly related to the order of the macroscopic gradients 
in the system. This was first pointed out by Grad in his monumental work on 
this subject [19] . Referring to the quantity that appears in his case namely, 
the one resulting from the collision integral in the thirteen moment solution to 
the Boltzmann equation, he asserts that "such a parameter has a number of 
interpretations" (see p. 271-72 of Ref. [19]) . Further, if one wants to classify 
the resulting equations in terms of power in the gradients one must apply the 
Chapman-Enskog expansion not to the Boltzmann equation but to the moment 
equations. It is only then, as has been extensively discussed in Ref. [20] one 
may extract the Euler, Navier-Stokes-Fourier, Burnett and higher order in the 
gradients contributions. Therefore Balescu's nearly exact results as he claims, 
obtained with the 21 moment method, are not clearly related to the ordinary 
hydrodynamic hierarchy of equations. We believe this is the main reason why 
our results are somewhat different from his. It is outside the scope of this paper 
to attempt a detailed comparison of both methods, mainly a purely algebraic 
issue. 

Finally we wish to remark also that these results together with similar ones 
here obtained for the Dufour coefficient may be useful in accounting for dissi- 
pative phenomena which are becoming rather important in the physics of the 
intracluster medium [21J-|23J. 
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1 2 3 4 5 6 

Fig. 1. The three thermal conductivities as functions of x — ui e T for T = 10 6 K 
and n = lO 15 ™" 3 . 
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